The effect of noise on the dynamics of a complex map at the 

period-tripling accumulation point 



Olga B. Isaeva, Sergey P. Kuznetsov 
and 

Tt - ■ Andrew H. Osbaldestin 

O 
O 
<N 

bJ[)! Institute of Radio -Engineering and Electronics of RAS, Saratov Branch, 



=5 
< 



u 

d 



Zelenaya 38, Saratov, 410019, Russia 



Department of Mathematics, University of Portsmouth, 
Portsmouth, P01 3HE, UK 



Abstract 



As shown recently (O.B. Isaeva et al., Phys.Rev E64, 055201), the phenomena intrinsic to 
dynamics of complex analytic maps under appropriate conditions may occur in physical sys- 
tems. We study scaling regularities associated with the effect of additive noise upon the period- 
tripling bifurcation cascade generalizing the renormalization group approach of Crutchfield et al. 
(Phys.Rev. Lett., 46, 933) and Shraiman et al. (Phys.Rev. Lett., 46, 935), originally developed 
for the period doubling transition to chaos in the presence of noise. The universal constant 
determining the rescaling rule for the intensity of the noise in period-tripling is found to be 
' 7 = 12.2066409... Numerical evidence of the expected scaling is demonstrated. 

PACS numbers: 05.45.-a, 05.45.Df, 05.10.Cc, 05.40.Ca 

1 INTRODUCTION 

p< Application of the renormalization group (RG) analysis and the concepts of universality and scaling in 
nonlinear dynamics started with the works of Feigenbaum concerning the period-doubling transition 
^ . to chaos pp. The universal quantitative regularities intrinsic to this type of behavior are common 
for a wide class of systems including one-dimensional maps, forced dissipative nonlinear oscillators, 
Rossler and Lorenz models, experimental systems in hydrodynamics, electronics, laser physics, etc 
In the context of real systems, a question of vital importance is the effect of noise on the dynamics 
at the onset of chaos. It is more or less obvious that the presence of noise destroys subtle details of 
the small-scale (or long-time) dynamics. The RG analysis developed by Crutchfield et al. [3] and by 
Shraiman et al. jl] gives a quantitative measure for this effect. Namely, to observe one more level 
of period doubling it is needed to decrease the noise amplitude by the factor 7 = 6.61903. 

After Feigenbaum the RG approach was developed by many authors, e.g. for the onset of chaos 
via quasiperiodicity and intermittency as well as for some situations arising in the multipa- 

rameter analysis of transition to chaos [7j or to strange nonchaotic attractors jHJ. In particular, 
the effects of noise on the dynamics have been studied in dissipative systems for intermittency [Oj, 
quasiperiodicity ^0], bicritical behavior jTT], and in Hamiltonian systems for period doubling ^2] 
and KAM-torus destruction [TT?] . 
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A special kind of dynamics occurs in iterative complex analytic maps z n+ i = f(z n ) [14J. An 
example is a complex quadratic map 

z n +i = A — z 2 n , (1) 

which gives rise to remarkable fractal formations in the complex plane of the variable z (the Julia 
sets) and of the complex parameter A (the Mandelbrot set). The latter is defined by the condition 
that being launched from the origin, the iterations at a given A never diverge. The Mandelbrot set 
has a subtle and complicated structure, which is a subject of extensive research. In particular, it 
contains special points of accumulation of bifurcation cascades. Beside the period doubling, there 
are also cascades of period tripling, quadrupling, etc. The period-tripling accumulation point has 
been studied first by Golberg, Sinai, and Khanin [THj (see also [IB]), and will be referred to as the 
GSK critical point . In the map (JTJ) it is located at 

A c = 0.0236411685377 + 0.78366065080521 (2) 

By analogy with Feigenbaum universality, it could be expected to occur in other nonlinear systems 
as well. 

However, in fact, complex analytic functions represent a very special and restricted class of maps 
because real and imaginary parts of f(z) must satisfy the Cauchy-Riemann equations. If this is not 
the case, the dynamics become drastically different In particular, it was shown that the period- 
tripling type of behavior does not survive, in general, under a non-analytic perturbation (in 
spite of the claim in the original paper [T3] and some other works jHH^Hl)- So, a principal question 
is: Do the phenomena intrinsic to complex analytic maps have any concern to dynamical behavior 
of physical systems? One example of an appropriate physical system was suggested by Beck fl\\. 
Another approach developed in Ref. [22] is based on a construction of a system of two coupled units, 
each of which can demonstrate Feigenbaum's period-doubling cascade, and the coupling is of some 
special kind. Moreover, this idea has been verified in an experiment with coupled electronic circuits. 

Recognizing a possibility of the physical occurrence of the phenomena of complex analytic dy- 
namics, we intend to study in this paper the effect of noise on the period-tripling cascade in the spirit 
of the earlier works of Crutchfield et al. [Sj and Shraiman et al. jl], on a basis of the appropriate 
generalization of the RG approach. 

2 RENORMALIZATION GROUP ANALYSIS 

If we perform 3 k iterations of the map at the GSK critical point A = A c , the renormalized evolution 
operator will converge, as known, to a universal function 

\imf k (zf k (0))/f k (0) = g(z), (3) 

fe^oo 

which satisfies the fixed-point RG equation 

g(z) = ag(g(g(z/a))). (4) 

Obviously, g(z) represents an evolution operator at the critical point for an asymptotically large 
number of iterations 3 k in terms of the properly normalized dynamical variable. Numerically, the 
approximation for this function was found in Ref. ^H] as a finite Taylor series: 

g( z ) 1 + (0.0547 - 0.7490z)^ 2 + (-0.0244 - 0.0525z)z 4 + ... (5) 
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together with the complex scaling constant 

a = l/g(g(g(Q))) ^ -2.09691989 + 2.35827964i (6) 
Let us introduce noise and consider a stochastic equation 

z n+1 = g(z n ) + eip(z n )£ n , (7) 

where e is a small parameter of the noise intensity, (p(z) is a smooth real function of complex argument, 
and £ n is a complex stationary random sequence with statistically independent subsequent terms. 
We assume that it has zero mean, unit mean square = 1, and zero correlation of the real 

and imaginary parts, ((Re£ n )(Im£ n )) = 0. It is expected that the scaling properties of the response 
under study will be independent on the concrete form of the distribution function, although it is 
convenient to suppose that a distribution function for £ n is of such kind that the amplitude of the 
noise is bounded. 

Three-fold application of the stochastic map yields, in the first order in e, 

z n +3 = g(g(g(zn))) + £{g'(g(g(zn)))g'(g(zn))<p(z n )£n , g x 
+g'(g(g(zn)))<p(g(zn))€n+i + <p(g(g(z n )))£ n+2 }. 

Now, we renormalize the variable z by substitution z — > z/a, and obtain 

z n +3 = ag(g(g(z n /a))) + ea{g'(g(g(z n /a)))g'(g(z n /a))ip(z n /a)£ n 
+g'(g(g(zn/a)))ip(g(z n /a))£ n+1 + <p(g(g(z n /a)))£ n+2 }. 

The first term in the right-hand part of the equation equals g(z), in accordance with Eq. Con- 
cerning the remaining terms, we make the following important remark. Let us suppose that we start 
at some z n . Consider an ensemble of the random numbers {£ n , £ n +i, Cn+2} an d compose them with 
complex coefficients given by functions of z n . As {Cm Cn+i ; Cn+2} are independent, the sum can be 
represented again as a random complex number with zero mean and unit mean square multiplied by 
a real function of complex argument, namely, 

a{g'(g(g(z n /a)))g'(g(z n /a))ip(z n /a)^ n ^ 
+g'{g(g(zn/a)))(p{g(z n /a))£ n+1 + <p{g{g{z n /a)))£ n+2 } = fi(z n )£ n . 

So, we rewrite © in the form analogous to (JJJ), with redefined function and random variable in the 
stochastic term: 

Zn+3 = g{z n ) + e<f(Zn)£ n . (11) 

To obtain a closed functional equation, we multiply both parts of Eq. (fTUj) by the complex conjugates, 
and perform an averaging over an ensemble of realizations of the noise. As = = 1, and 

(£n£n+i) = = O5 we come to the relation 

H*n) = W?{\g'{g{g(zn/oc)))g'{g{z n /a))\ 2 ${z n /a) + \g'{g{g{z n /u)))\ 2 <&(g(z n /a)) + ®(g(g(z n /a)))}, 

(12) 

where $(2) = [v?(^)] 2 . Obviously, Eq. (|T2|l has a structure $(z) = L$(z), where L is a linear operator 
of the functional transformation given by the right-hand part of (|12|) . Repetitive application of the 
same procedure to Eq. (|TT|) yields a sequence of functions with asymptotic behavior &k(z) — fl k Q(z), 
determined by eigenvector $(-2), associated with the largest eigenvalue Q of the operator L: 

£l$(z) = \a\ 2 {\g'{g{g{z/a)))g'{g{z/a))\^{z/a) + \g'(g(g(z/a)))\^(g(z/a)) + <S>(g(g(z/a)))}. 

(13) 
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Figure 1: Plot of the eigenfunction responsible for the effect of noise at the period tripling accumu- 
lation point obtained numerically (see text). 

As mentioned above, the universal function g(z) has been approximated as a finite expansion over 
even powers of the argument j!5j . Using this data, we have numerically constructed the functional 
transformation of the right-hand part of Eq ()13j) . The unknown function $(z) is represented by a set 
of real values at nodes of a grid in a square { — 1.2 < Kez, Im,2 < 1.2}, and by an interpolation scheme 
between them. Taking random initial conditions for $(-2), we perform the functional transformation 
and normalize the resulting function as = &(z)/&(0). This operation is repeated many times, 

until the form of the function stabilizes (see Fig. 1). The value of $(0). (before the normalization) 
converges to the eigenvalue Q = 149.0020828. 

In the linear approximation with respect to the noise amplitude, the stochastic map for the 
evolution over 3 fc steps at the GSK critical point asymptotically may be written as 

Z n +3" = g{Zn) + £J k( p(Zn)L, (14) 

where 

cp(z) = yj$°(z), 7 = = 12.20664093. (15) 

If we consider a small shift of A from the GSK critical point, then an additional perturbation 
term appears in the equation: 

z n+3 k = g(z n ) + Ci6ihi(z n ) + ej k <p(z n )£ n . (16) 

Here, h\{z) represents an eigenvector of the linearized RG equation without noise associated with 
the eigenvalue 5\ = 4.60022558 — 8.98122473z. The coefficient C\ depends on the parameter and 
vanishes at the critical point GSK. In a close neighborhood of the critical point it is sufficient to 
consider only the leading terms of the expansion and set C\ oc (A — A c ). 

Now, we are ready to formulate the basic scaling property that follows from (|16|) . 

If we triple the number of time steps (i.e., change k to k + 1), decrease the parameter difference 
A A = A — A c by division by Si, and decrease noise amplitude e by factor 7, then the form of the 
stochastic map (fTljj) remains unchanged. Thus, with the new parameters, (AA/5i, e/7) the noisy 
system will demonstrate the same behavior as that with the old ones, but with a tripled time scale. 
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3 MODEL MAP AND NUMERICAL EXPERIMENTS 



To verify the RG results in numerical experiments, we use a model map with additive noise: 

z n +i = A - zl + s£ n . (17) 

Real and imaginary parts for each term of the complex random sequence £ n , are obtained as sums of 
N = 10 zero-mean computer generated pseudo-random numbers, which are properly normalized to 
have (\C,n\) = 1- (I n accordance with the central limit theorem, £ n is very close to a Gaussian random 
variable.) 

Fig. 2 illustrates fine structure of the Mandelbrot set near the period-tripling accumulation point 
as it looks (a) in the absence and (b) in the presence of noise. The pictures represent the Lyapunov 
charts for the complex parameter A. (See |23| ITH] for previous applications of this method.) At each 
pixel of the two-dimensional plot we estimate the Lyapunov exponent 

N 

A = lim N' 1 y2 In \ 2z n \ (18) 

n=l 

from numerical computations, and mark the pixel with a gray tone. The tones vary from dark to light 
as the Lyapunov exponent varies from large negative to zero: white corresponds to zero, and black 
to positive Lyapunov exponent. Divergence is shown by uniform coloring with one special gray tone. 
The GSK critical point is located exactly at the center of the diagram. A small box containing the 
critical point, is magnified and rotated in accordance with multiplication by 8\ ~ 4.6002 — 8.9812*. 
As the parameter rescaling is associated with a tripling of characteristic time scale, we accompany 
it with reducing an interval for gray coding of Lyapunov exponents by factor 1/3 for each successive 
picture in the row. On the panel (a) such diagrams are shown for zero amplitude of noise, e = 0, 
and on panel (b) for noise of fixed intensity, e = 0.001. Observe the similarity of the leaves of the 
"Mandelbrot cactus" in the first case, and the progressive destruction of the structure at subsequent 
levels of the resolution in the second. 

In Fig. 3 we show analogous pictures with noise, but now its intensity is reduced for each subse- 
quent picture in a row by the factor 7 = 12.2066... found from the RG analysis. Panels (a) and (b) 
correspond to different initial levels of the noise, respectively, e = 0.003 and 0.009. Now, the simi- 
larity of the pictures is restored. With noise reduction by 7, just one more level of smaller leaves of 
the Mandelbrot cactus reveals itself near the GSK critical point. 

Fig. 4 illustrates similarity in structure of the noisy attractors for the model map (|17|) in the plane 
of the complex variable z. Diagrams (b) and (c) relate to noise intensities reduced by the factors 
7 and 7 2 in comparison with the plot (a). The parameter A values correspond to the superstable 
cycles of period 3, 9, and 27, respectively. The right-hand panels show the indicated parts of the 
pictures drawn with rescaling by complex factors a (diagram (b)) and a 2 (diagram (c)). Observe 
the similarity of the right-hand diagrams. 

4 CONCLUSION 

In this paper we have studied the effect of noise on a complex analytic map at the period-tripling 
accumulation point. It was shown that the effect of noise obeys regularities of a similar nature to 
those reported for period doubling in real one-dimensional maps. Namely, to observe one more level 
of the bifurcation cascade it is needed to decrease the intensity of noise by some definite scaling 
factor. For the case of the period-tripling this factor has been found to be 12.2066..., as follows from 
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Figure 2: The Lyapunov charts for the model map (|T7jl on the complex parameter plane for zero 
intensity of noise e = (a) and in presence of noise, e = 0.001 (b). The GSK critical point is 
located exactly at the center of the diagrams. The tones vary from dark to light as the Lyapunov 
exponent varies from large negative to zero: white corresponds to zero, and black to positive values. 
Divergence is shown by uniform coloring with one special gray tone. A small box containing the 
critical point, is magnified and rotated in accordance with multiplication by 5\ ~ 4.6002 — 8.9812z. 
The interval for coded values of Lyapunov exponents is reduced with factor 1/3 for each successive 
picture in the row to visualize self-similarity of the pictures in absence of noise. 
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Figure 3: Lyapunov charts for the stochastic model (|17|) on the complex parameter plane in presence 
of noise for e — E\ — 0.003, e = £1/7, and e = ei/7 2 (a), and for e = e<i = 0.009, e = £2/1, 
and e = E2/1 2 (b). The GSK critical point is located exactly at the center of the diagrams. The 
coordinate rescaling rule and the gray coding are the same as in Fig. 2. The interval for coded values 
of Lyapunov exponents is reduced with factor 1/3 for each succesive picture in the row to visualize 
the similarity of the pictures. g 
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Figure 4: Similarity in structure of the noisy attractors for the model map (jl7j) on the plane of 
complex variable z. The parameter A values are 0.122561167 + 0.744861767z (a), 0.031552975 + 
0.790783175* (b), 0.023369416 + 0.784679678* (c), and the noise intensities are e = 0.015, 0.015/7, 
and O.OI5/7 2 , respectively, where 7 = 12.2066 The right-hand pictures represent the definite parts 
of the attractors drawn with rescaling by complex factors a = —2.0959... + z2.3582... (b) and a 2 (c). 
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the renormalization group analysis. Also we have demonstrated the scaling properties associated 
with the new constant in numerical experiments. 

The undertaken analysis is essential for discussion of the possibility of observation of phenomena 
of complex analytic dynamics in real physical experiment, e.g. in mechanics or electronics [2UE21EB]- 
The presence of noise in such systems is inevitable, and our result gives a quantitative foundation 
for estimates of an observable number of levels in the bifurcation cascade. 

Our work puts the phenomena in complex analytic map into the same context in respect to the 
effect of noise, as other situations of universal scaling behavior: period doubling [3J Hj, intermit- 
tency quasiperiodicity (THj, bicriticality [TTj, scaling phenomena in Hamiltonian systems [T2| IT3]. 
So, this specific field is enriched with one more nontrivial example of the scaling behavior linked with 
presence of noise. 

The approach developed in this article may be adapted to study the effect of noise on other 
bifurcation cascades in complex analytic maps ("n-tupling" with n = 4,5, ... fHJE3)- Moreover, it 
gives a basis from which to pose more general questions on the effect of noise on the Mandelbrot-like 
sets in physical systems possessing such kind of objects in their parameter space. 
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